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INTRODUCTION 


Propagation  of  ground-based,  high-energy  laser  (HEL)  beams  through  the  atmosphere 
at  a  laser  wavelength  of  A l  =  1.06  //m  is  limited  by  thermal  blooming  due  to  aerosol 
absorption.  The  scaling  of  critical  laser  power  with  beam  diameter  requires  that  the 
aperture  size  exceed  20  m  for  average  laser  powers  of  100  MW. 

Science  Research  Laboratory  (SRL)  has  considered  the  resolution  of  this  problem  by 
vaporizing  the  aerosols  with  a  designated  laser  for  aerosol  vaporization  (LAV).  In  this 
report,  the  dependence  of  vaporization  on  fluence  and  wavelength  will  be  described.  A 
laser  wavelength  A„  =  1.06  //m  and  intensity  of  10  MW/cm2  yields  a  reduction  in  the 
atmospheric  absorption  by  a  factor  of  5  at  2  //sec,  corresponding  to  20  J/cm2  of  laser 
fluence.  The  best  results  are  obtained  at  A„  =  0.248  // m,  which  is  the  KrF  laser  wavelength. 
In  this  case,  an  intensity  of  3  MW/cm2  and  a  pulse  duration  of  1  //sec,  corresponding  to  a 
fluence  of  3  J/cm2,  yields  a  reduction  by  a  factor  of  50.  Such  a  decrease  in  the  absorption 
coefficient  eliminates  the  difficulties  associated  with  thermal  blooming. 

The  dependence  of  the  critical  HEL  power  on  the  beam  aperture  diameter  d  has  been 
determined  in  the  presence  and  absence  of  LAV,  assuming  no  wind  shear.  Without  vapor¬ 
ization  the  scaling  law  is  pessimistic — one  can  propagate  less  than  10  MW  average  power 
at  d  =  10  m.  If  one  uses  a  KrF  laser  at  an  average  power  of  2  MW,  a  pulsewidth  of  1 
//sec,  and  a  pulse  repetition  frequency  (PRF)  of  10  Hz,  the  critical  power  at  d  =  5  m  is 
100  MW.  This  greatly  alleviates  the  thermal-blooming  problem. 
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1.0  AMELIORATION  OF  THE  THERMAL  BLOOMING  PROBLEM 

The  main  obstacle  for  beam  propagation  of  ground  based,  high  energy  lasers  (HEL)  is 
the  thermal  blooming  in  the  atmosphere  due  to  molecular  and  aerosol  absorption.  Even 
a  small  fraction  of  the  laser  energy  deposited  in  the  air,  induces  changes  in  the  refractive 
index  resulting  in  a  significant  beam  spreading.  This  phenomena  is  referred  to  as  thermal 
blooming.  The  laser  wavelength  is  chosen  in  an  atmospheric  window  where  the  absorption 
is  small.  For  BMD  application  the  current  choice  is  A  =  1.06  /zm.  The  laser  can  be 
tuned  to  a  wavelength  where  there  are  no  strong  water  vapor  absorption  lines  and  aerosol 
absorption  dominates. 

The  HEL  beam  path  in  the  atmosphere  can  be  cleared  by  a  designated  laser  for  aerosol 
vaporization  (LAV).  Such  a  laser  operates  at  an  average  power  only  a  few  percent  of  the 
HEL  power.  To  minimize  the  LAV  power  one  has  to  choose  a  wavelength  at  which  the 
aerosol  absorption  coefficient  (a)  is  large.  At  A  =  1  /z m,  a  a  10-7  cm-1,  while  at  A  = 
0.248  fj,m,  a  ~  2  x  10-6  cm-1.  As  a  result  the  vaporization  fluence  for  a  KrF  laser  is 
nearly  20  times  smaller  than  for  a  Nd:glass  laser.  To  vaporize  the  aerosols  at  A  =  0.248 
fxm  one  needs  a  laser  fluence  >  1  J/cm2  and  in  this  regime  aerosol  absorption  dominates 
the  saturated  ozone  absorption  in  the  lower  atmosphere. 

The  temporal  format  of  the  LAV  should  be  consistent  with  the  heat  exchange  time 
between  the  aerosols  and  the  surrounding  air,  as  well  as  the  wind  velocity.  For  typical  dry 
aerosols  the  air  is  heated  at  times  >  1  ftsec  and  the  LAV  should  vaporize  the  particles  in  less 
than  1  fisec.  This  determines  the  LAV  pulse  length.  If  the  vaporization  laser  propagates 
in  a  guard- ring  around  the  HEL  beam,  its  pulse  repetition  frequency  (PRF)  is  determined 
by  the  thickness  of  the  guard-ring  and  the  wind  velocity.  Typical  parameters  lead  to  a 
PRF  of  10  Hz. 

To  reduce  the  thermal  blooming  and  propagate  100  MW  of  HEL  power  at  beam  diam¬ 
eters  in  the  range  of  4m,  one  needs  only  2  MW  of  LAV  power.  This  is  an  excellent  trade-off 
if  one  considers  that  in  the  absence  of  LAV  only  10  MW  of  power  can  be  transmitted  to  a 
relay  mirror  in  space.  The  alternative  to  LAV  is  either  a  very  large  and  expensive  beam 
director  ( d  >  20  m)  or  a  risky  and  unproven  technology  based  on  optical  phase  conjugation, 
i.e.  Brillouin  scattering  or  four-wave  mixing. 

Science  Research  Laboratory  will  propose  to  demonstrate  the  clearing  of  aerosols  by  a 
laser  for  aerosol  vaporization  in  a  series  of  laboratory  experiments: 

•  Single  pulse  LAV 

-  Determine  the  changes  in  the  aerosol  size  distribution  and  the  reduction  in  aerosol 
mass  loading. 

-  Determine  the  reduction  in  the  aerosol  absorption  coefficient  after  vaporization. 

•  Repetitively  pulsed  LAV 

-  Determine  the  reduction  in  thermal  blooming  in  steady  state. 
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2.0  THEORY  OF  AEROSOL  VAPORIZATION 
2.1  INTRODUCTION 

Propagation  of  high  power  laser  beams  through  the  atmosphere  gives  rise  to  molecular 
and  aerosol  absorption  which  results  in  local  heating.  The  air  expands  to  restore  pressure 
balance  and  the  resulting  density  perturbation  leads  to  a  decrease  in  the  local  refractive 
index.  The  phase  front  of  the  HEL  beam  moves  faster  in  regions  of  lower  refractive  index 
(n)  and  bends  toward  regions  of  larger  n.  The  result  is  a  thermal  blooming  of  the  beam  in 
the  cross  wind  direction  and  a  tilt  in  the  wind  direction/1-2)  The  propagation  of  a  ground 
based  laser  beam  to  a  relay  mirror  in  space  requires  compensation  for  the  deleterious 
effects  of  thermal  blooming  and  atmospheric  turbulence.  This  is  achieved  with  a  closed 
loop  adaptive  optics  system/3-5)  A  beacon  in  space  is  used  to  probe  the  atmosphere  and 
sense  index  of  refraction  perturbations  which  are  measured  at  the  transmitter  on  the 
ground.  Phase  compensation  means  that  the  conjugate  of  the  incoming  beacon  phase 
be  imposed  on  the  outgoing  laser  beam  with  a  deformable  mirror  consisting  of  actuators 
spaced  a  coherence  length  r0  apart.  Such  a  system  compensates  very  well  for  atmospheric 
turbulence  with  Strehl  ratios  in  excess  of  80%.  The  thermal  blooming  can  be  compensated 
only  up  to  a  critical  laser  power,  Pc.  The  scaling  law  of  Pc  as  a  function  of  beam  diameter 
(d),  laser  wavelength  (A)  and  absorption  coefficient  (ck)  is  of  the  form,  Pc  ~  da X^/a.  The 
values  of  a  and  b  are  the  subject  of  an  active  study  with  no  consensus  at  the  present 
time.  The  fact  that  Pc  is  inversely  proportional  to  a  is  beyond  dispute.  It  is  clear,  that 
despite  uncertainties  in  the  scaling  law,  the  propagation  of  100  MW  of  average  GBL  power 
necessary  for  BMD  applications  will  result  in  unacceptably  large  aperture  sizes.  The  only 
effective  remedy  is  to  reduce  significantly  the  absorption  coefficient  a. 

The  laser  wavelength  is  chosen  in  an  atmospheric  window,  where  a  due  to  molecular 
and  aerosol  absorption  is  minimum.  One  of  the  atmospheric  windows  is  near  1.06  /j, m. 
At  this  wavelength  the  molecular  absorption  is  dominated  by  the  water  vapor  continuum. 
Experiments  at  SRL  in  the  past  few  months  have  demonstrated  that  for  A  ~  1.057  fim 
there  is  a  micro-window  of  10  cm-1  width  where  a  is  of  the  order  of  10-9  cm-1.  The  results 
are  shown  in  Figs,  la  and  lb.  In  Fig.  lb,  the  continuous  line  is  a  fit  from  FASCODE 
and  the  dots  are  the  SRL  data.  By  comparison,  a  from  aerosol  absorption  in  the  1  fim 
range  and  good  visibility  (23  km)  is  of  the  order  of  10-'  cm-1.  If  the  GBL  wavelength 
is  in  the  water  vapor  absorption  window,  the  aerosols  are  by  far  the  dominant  absorber. 
To  reduce  a  we  consider  clearing  the  aerosols  by  vaporization  with  a  designated  laser  for 
aerosol  vaporization  (LAV). 

The  particles,  responsible  for  the  absorption,  range  in  size  between  0.1  fj, m  and 
0.5  /mi/6)  As  they  are  heated  by  the  laser  they  transfer  the  heat  to  the  surrounding 
air  on  a  time  scale  of  1  /t sec  -  10  /tsec/7-8)  Efficient  vaporization  takes  place  when  the 
laser  heats  the  particles  to  the  boiling  point  in  a  short  time  of  the  order  of  1  fi sec  and 
this  consideration  determines  the  pulselength  of  LAV.  The  pulse  format  of  the  induction 
linac  HEL  is  such  that  a  single  50  nsec  pulse  can  not  heat  the  aerosols  to  the  vaporization 
temperature.  With  a  PRF  of  the  order  of  10  kHz  the  particles  cool  down  between  pulses 
and  do  not  accumulate  the  heat.  This  is  the  main  reason  why  the  IL-HEL  cannot  vaporize 
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760  TORR,  22.9°C 
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Figure  la:  Atmospheric  absorption  near  1.06  pm  (FASCODE) 


ABSORPTION  COEFFICIENT  (1/cm) 


Figure  lb:  Water  vapor  absorption  in  N2 .  T=30°C,  vapor  pressure  14  Torr, 

total  pressure  767  Torr.  The  discrete  points  are  SRL  measurements. 
The  solid  line  is  a  FASCODE  result  convolved  with  the  SRL 
line  shape  and  normalized  to  a  known  line  strength. 
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the  aerosols.  The  heating  of  the  particles  is  related  also  to  the  imaginary  part  of  their 
refractive  index,  n,\  n,-  is  a  sensitive  function  of  A.  The  HEL  wavelength  is  chosen  to 
minimize  n,-  and  the  absorption,  while  that  of  the  LAV  should  be  chosen  at  A„  where  n,-  is 
large.  From  our  calculations  it  appears  that  \v  =  0.248  /tm  is  the  optimum  wavelength. 
This  will  insure  that  the  power  of  the  LAV  beam  is  a  small  fraction  of  the  power  of  the 
BMD-HEL  beam.  If  one  considers  a  LAV  beam  propagating  in  a  guard-ring  of  thickness 
dv  surrounding  the  HEL  beam,  the  PRF  will  be  determined  by  Vw/dv,  where  Vw  is  the 
wind  velocity.  With  Vw  =  5  m/sec  and  dv  =  0.5  m  one  finds  PRF  =  10  Hz. 

In  this  section  the  dynamics  of  single  particle  heating  and  vaporization  is  studied,  and 
the  critical  laser  fluence  is  determined.  The  change  in  size  distribution  and  absorption 
coefficient  during  aerosol  vaporization  are  found  as  function  of  fluence.  The  turbulence 
induced  intensity  scintillations  of  the  LAV  lead  to  phase  distortions  on  a  scale  of  the  order 
of  the  coherence  length  of  the  LAV.  Their  effect  on  the  propagation  of  the  HEL  beam  is 
studied.  The  basic  requirements  for  LAV  are  determined  and  its  impact  on  the  thermal 
blooming  scaling  and  the  overall  system  are  evaluated. 

2.2  HEATING  AND  VAPORIZATION  OF  THE  AEROSOLS 

The  atmospheric  aerosols  have  components  of  different  chemical  composition,  index  of 
refraction  and  sizes.  The  large  particles  with  dimensions  above  1  p m  consist  mostly  of  sand 
and  are  irregular  in  shape.  However,  their  number  is  quite  small  (<  100  particles/cm3) 
and  their  imaginary  index  of  refraction,  n,-,  is  small.  Their  contribution  to  the  absorption 
coefficient  is  considered  negligible.  The  very  small  particles  with  a  <  0.05  pm  are  probably 
mostly  carbon  with  large  nj,  but  their  mass  loading  is  negligible  and  their  contribution 
to  a  will  be  considered  small.  The  particles  which  are  responsible  for  thermal  blooming 
typically  are  in  the  range  between  0.05  jum  and  0.5  jun.  We  shall  assume  that  they  are 
uniform  and  spherical  in  shape.  Since  the  mean  free  path  in  the  atmosphere  is  0.05  fj,m,  one 
treats  the  process  of  heat  transfer  from  the  particles  to  the  air  as  a  conduction  dominated 
diffusion  process,  i.e.  the  temperature  in  the  air  surrounding  the  particles  is  given  by  the 
heat  diffusion  equation/7-8) 


pa,Ca,Ka  are  the  density,  specific  heat  and  thermal  conductivity  of  the  air  and  in  our 
model  they  are  assumed  to  be  constant.  The  temperature  inside  the  particle  is  considered 
uniform  because  its  thermal  conductivity  is  much  greater  than  that  of  the  air.  Typical 
temperature  equilibration  times  for  particles  of  1  pm  radius  are  of  the  order  of  10  nsec. 
Furthermore,  the  laser  intensity  inside  the  aerosol  is  assumed  to  be  uniform  and  the  effect 
of  focusing  of  the  light  inside  the  particle  is  neglected.  The  conditions  of  a  diffusion 
process  and  no  focussing  lead  to  larger  threshold  values  for  vaporization.  Anomalously  low 
intensity  threshold  values  have  been  observed  for  particles  with  real  part  of  the  refractive 
index  approaching  2,  where  focusing  inside  a  sphere  is  occurring.  In  general,  our  analysis 
will  conservatively  estimate  the  vaporization  threshold. 
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The  boundary  condition  at  the  particle  surface  is 

T  ^7r  3/--T  &T  2  Tf 

I°A  =  -Ja  Crpr~di  ~  ^  Kad 7 


4:‘KOiLvpv 


da 

dt 


(2) 


a  a  is  the  aerosol  absorption  cross  section  calculated  form  Mie  theory/9)  pp,Cp,Lv  are  the 
particle  density,  specific  heat  and  latent  heat  of  vaporization,  and  I  is  the  laser  intensity.  To 
determine  the  vaporization  threshold,  consider  T  just  below  the  boiling  point  and  neglect 
the  last  term  in  Eq.  (2).  Equation  (1)  with  the  condition  (2)  was  solved  analytically  for  a 
CW  laser/8)  The  result  can  be  extended  to  a  laser  pulse  length  rp.  At  the  particle  surface 
one  finds 


T(a,  t)  =  T0  + 


I*A 

AvaK, 


m 


1  —  exp 


(-1)]  -«.(*-  rr)  [l  -  exp  (-^)]  }  .  (3) 


where  T0  is  the  ambient  temperature,  ra  =  a2ppCp/3Ka  is  the  equilibration  time  and  8(t ) 
is  the  unit  step  function,  8{t)  —  1  when  t  >  0  and  zero  otherwise.  The  absorption  cross 
section,  a  a  is  a  function  of  aerosol  radius,  a,  the  laser  wavelength,  A  and  the  index  of 
refraction,  n  =  nr  +  ini .  The  latter  is  shown  as  a  function  of  A  from  0.248  pm  to  10 
pm  for  representative  dustlike  and  water-soluble  aerosols.  Note  from  Fig.  2  that  the  real 
part  is  between  1  and  1.5  for  all  wavelengths,  except  the  range  9-10  pm,  where  focussing 
inside  the  aerosol  is  possible.  The  imaginary  part,  see  Fig.  3,  has  peaks  near  the  shortest 
and  longest  wavelengths.  The  absorption  efficiency,  Qa  =  &  a/ 7T  a? ,  is  calculated  from  Mie 
theory  applied  to  dry  aerosols  at  A  =  0.248  pm,  1.06  pm.  The  result  is  shown  on  Fig. 
4.  The  regime  at  small  radii,  where  Qa  oc  a  and  a  a  oc  a 3  is  called  the  Rayleigh  regime 
and  for  calculating  absorption  it  is  determined  by  a  <  8,  where  8  =  (2fcra;)-1  is  the  bulk 
absorption  distance.  At  A  =  1.06  pm,  n,-  =  8  x  10— 3,  one  has  8  =  10  pm  and  the  Mie 
theory  result  falls  into  the  Rayleigh  regime  for  a  <  5  pm.  At  A  =  0.248  pm,  n,-  =  3  x  10— 2 , 
8  =  0.7  pm  and  the  Rayleigh  regime  applies  to  particles  of  sizes  up  to  0.4  pm.  Since  the 
submicron  particles  dominate  the  absorption  we  treat  the  aerosols  in  the  Rayleigh  regime 
and  write  for  a  a, 

7rC[3  fA\ 

°A  =  —  .  (4) 

The  vaporization  threshold  is  evaluated  from  Eq.  (3)  when  T(a,  rp)  —  Tv,  the  vaporization 
temperature.  In  Fig.  5  the  critical  fluence,  Fc  —  Icrp,  is  plotted  as  a  function  of  tp  for  Tv 
=  2500°C,  different  laser  wavelengths  and  dry  aerosols  of  the  dominant  size  a  =  0.2  pm. 
When  Tp  <  r0  the  critical  fluence  is  independent  of  pulselength.  From  Eq.  (3)  one  can 
write  ( tp  «  r0), 

r„  —  Tp  47m3 

rc—  o  Pp'-'P  * 

C  a  & 

with  a  a  in  the  Rayleigh  regime  and  Tv»  T0  one  finds 


Fr  — 


\pp  C p'Fy 

37rn; 


(5) 


(6) 


As  seen  from  Fig.  5,  Fc  is  independent  of  rp  for  rp  <  lpsec  because  the  equilibration 
time  for  0.2  pm  aerosols  is  ra  =  1.3  psec.  The  threshold  is  proportional  to  A/n/A)  and 
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this  quantity  determines  the  effectiveness  of  LAV.  Clearly,  the  KrF  wavelength  is  the  most 
effective  with  Fc  =  0.5  J/cm2.  In  Fig.  6  the  critical  fiuence  Fc(tp)  is  plotted  for  dry 
aerosols  of  different  sizes  and  A  =  0.248  ^m.  Note  that  a  fiuence  of  1  J/cm2  at  rp  =  1 
//sec  will  vaporize  all  particles  larger  than  0.2  p m  in  radius.  Figures  5  and  6  illustrate 
also  that  lasers  of  longer  pulselenghts  require  substantially  larger  fluencies  of  vaporization 
because  the  particles  exchange  heat  with  the  air  during  the  laser  pulse.  For  example,  at  A 
=  1.06  /zm,  tp  =  10  //sec,  a  =  0.2  /zm  the  critical  fiuence  exceeds  100  J/cm2,  see  Fig.  5. 
This  important  result  will  be  used  in  the  experiments  proposed  by  SRL  to  combine  a  short 
pulse  LAV  and  a  long  pulse  probe  laser  of  large  fiuence  to  achieve  the  required  sensitivity 
in  measuring  the  absorption  coefficient  a. 


When  F  >  Fc  a  rapid  evaporation  takes  place  and  the  last  term  on  the  r.h.s.  of  Eq. 
(2)  dominates.  With  a  a  given  by  Eq.  (4)  one  can  write^10) 


da  irriil 

dt  XppLv 


(7) 


The  particle  is  reduced  in  size  exponentially, 


a  =  a0exp 


where  Fv  =  \pPLv/Z’Krii  is  the  fiuence  vaporization  rate.  For  dry  aerosols  typical  values 
axe  pp  =  2.3  g/cm3  and  Lv  =  4kJ/g. 


Table  Is  Vaporization  Fiuence  for  Dustlike  Aerosols  at 
Different  Wavelengths. 


LAV  Wavelength 
(pm) 

Refractive  Index 
of  Aerosols 

Vaporization  Fiuence 
(J/cm2) 

0.248 

1.52  +  3xl0-2  i 

0.8 

1.06 

1.52  +  8xl0“3  i 

12 

10.6 

1.62  +  0.12  i 

8 

The  reduction  in  particle  size  means  reducing  the  mass  loading  of  the  aerosols  and 
reducing  a.  To  determine  that  effect  we  study  the  whole  ensemble  of  particles  of  different 
sizes. 


2.3  SIZE  DISTRIBUTION  OF  THE  AEROSOLS  WITH  VAPORIZATION 

The  size  distribution  of  the  aerosols  has  been  the  subject  of  extensive  research  in  the 
past  and  it  is  continuing  at  the  present  time.  It  appears  that  the  aerosols  consist  of  at  least 
2  or  3  components  with  different  distributions,  ranging  in  size  from  10-2  pm  to  a  few  pm. 
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■o -  WATER-SOLUBLE 

- -  DUSTLIKE 


5  10 


WAVELENGTH  ([Jin) 


Figure  2:  Real  part  of  refractive  index  of  aerosols 
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Figure  3:  Imaginary  part  of  refractive  index  of  aerosols 
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Q  (ABSORPTION  EFFICIENCY) 


Figure  4:  Absorption  efficiency  of  dustlike  spherical  particles  for  Nd 
glass  and  KrF  laser  wavelengths 
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FLUENCE  (J /  cm 


Figure  5:  Critical  fluence  for  aerosol  evaporation 
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It  is  assumed  that  the  submicron  fraction  is  the  most  important  for  the  optical  properties 
of  haze  with  the  Yunge  and  modified  gamma  distributions  being  favored  by  researchers  in 


this  field.  While  the  Yunge  distribution,  / 


l,  fits  well  the  data,  it  suffers  from  the 


drawback  that  without  cutoffs  it  leads  to  divergent  moments.  For  theoretical  studies  the 
modified  gamma  distribution  function  is  preferable.  It  is  given  by*10) 


,,  .  .r  4a2  (  2 a\ 


where  a0  is  the  most  probable  size  of  the  aerosols  and  N0  is  the  number  of  particles/cm3, 
Jo°°  f(a)da  —  N0.  Equation  (9)  has  been  used  in  the  study  of  cloud  clearance  as  a  size 
distribution  for  the  water  droplets.  Experimental  results  with  a  CO2  laser  have  shown  an 
excellent  agreement  with  theory.  The  size  distribution  satisfies  the  equation*10) 

f  +  |(a/)  =  0  (10) 

It  is  assumed  that  during  vaporization  particles  are  not  broken  up  into  fragments  and 
they  do  not  disappear.  The  aerosols  are  only  reduced  in  size  and  their  number  (N0)  is 
conserved.  From  Eqs.  (8)  and  (10)  one  can  write  for  /  in  terms  of  fluence  and  particle 
radius, 

dF~Wvd^^  =  ° 

The  solution  of  Eq.  (11)  with  an  initial  distribution  (9)  is  given  by 


4  n? 

f(a,F)  =  N0— exp 


2  a 

exp - exp 


At  the  present  time  there  is  no  direct  measurement  of  a  size  distribution  function  with 
vaporization.  SRL  proposes  to  do  an  experiment  to  measure  /  before  and  after  vaporization 
and  this  will  provide  a  crucial  link  to  determine  the  vaporization  mechanism  and  find  the 
absorption  coefficient  a  as  a  function  of  fluence. 

2.4  ABSORPTION  COEFFICIENT  WITH  AEROSOL  VAPORIZATION 

The  absorption  coefficient  for  the  aerosols  is  determined  from  the  absorption  cross- 
section  and  the  size  distribution, 


00 

a(F)  =  J  a  A(a)f(a,F)dc 


From  Eqs.  (4)  and  (12)  one  obtains 


a(F)  =  ot0exp  ( 


where 


oc0  =  N0 


307 r2n,a3 
A 
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is  the  absorption  coefficient  in  the  absence  of  vaporization.  Equation  (14)  is  a  central  result 
of  this  proposal  and  shows  why  a  fluence  above  the  vaporization  rate  significantly  decreases 
the  aerosol  absorption.  Note  that  in  the  Rayleigh  regime  the  absorption  coefficient  is 
proportional  to  the  mass  loading  (a a  ~  a3)-  With  vaporizaiton  the  aerosol  mass  simply 
decreases  as  shown  by  (14),  which  is  consistent  with  the  exponential  decrease  in  radius, 
see  Eq.  (8).  This  means  that  at  A  =  0.248  fi m  and  fluence  of  4  J/cm2  the  absorption 
coefficient  is  reduced  by  two  orders  of  magnitude,  see  Table  1.  In  that  case  the  effects 
of  thermal  blooming  would  be  negligible.  The  exponential  decrease  of  a  as  a  function  of 
fluence  has  been  observed  experimentally  in  the  case  of  water  droplet  evaporation  for  cloud 
clearing. 

An  indirect  confirmation  of  the  validity  of  our  size  distribution  and  absorption  cross- 
section  models  comes  by  comparing  a0  from  Eq.  (15)  with  the  results  in  the  Handbook  of 
Optics.  There  a  Yunge  distribution  and  Mie  scattering  theory  have  been  used.  At  23  km 
visibility  N0  =  2.8  xlO3  particles/cm3.  At  A  =  1.06  /xm  for  dry  aerosols  (30%),  n,  =  8  x 
10-3  and  for  soluble  aerosols  (70%),  n,-  =  1.7  x  10"2.  With  a0  =  0.1  A<m  one  finds  from 
Eq.  (15)  a  =  1.2  x  10-7  cm-1.  The  result  in  the  Handbook  of  Optics  is  a  =  1.1  x  10-7 
cm-1.  The  table  below  shows  an  excellent  agreement  at  different  wavelengths  where  data 
is  available. 

Table  2:  Comparison  of  Absorption  Coefficients  at  Different  Wavelengths 


A(jum) 

a(1.06)/o:(A) 
Theory,  Eq.  (15) 

a(1.06)/a(A) 
Handbook  of  Optics 

10.6 

1.68 

1.97 

3.392 

4.95 

4.86 

0.6328 

1.37 

1.18 

0.337 

0.82 

0.74 

0.248 

0.12 

Not  available 

There  is  a  clear  distinction  between  the  absorption  coefficients  of  aerosols  and  the  air. 
The  heat  is  transferred  from  the  aerosols  to  the  surrounding  air  with  a  time  delay  given  by 
the  equilibration  time.  The  aa  which  enters  the  thermal  blooming  analysis  is  a  measure 
of  the  rate  of  heat  deposition  in  the  air/8) 

paCa  dT 

<Xa~~7~~di 

where  T  is  the  temperature  of  the  air.  aaI  is  the  rate  of  heating  per  unit  volume  of  air.  To 
establish  the  relationship  between  oca  from  (16)  and  a  from  (13)  let  us  write  the  increase 
of  the  temperature  in  the  air  by  solving  Eqs.  (1)  and  (2), 


A  T(r,rj,<,a,)  = 


i(nW(«.)  B 

4  vKa  |  r  —  r,  | 


(17) 
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where 


B(r  -  ?i,  t )  = 


r  -  r,-  |  — a,-  Kg 

y/Wi  ’  PaCa  * 

The  temperature  increase  at  point  f  and  time  t  is  due  to  an  aerosol  of  radius  a,-  at  a  point 
f,-.  The  air  temperature  due  to  an  ensemble  of  particles  is 


(AT)  =  ^AT(f,ri!i,ai)  .  (18) 

i 

If  the  intensity  variation  on  a  scale  comparable  to  the  interparticle  distance  is  small  and 
they  are  many  particles  in  a  volume  d3,  where  d  is  the  laser  beam  diameter,  one  can  write 
(AT)  in  the  form  of  an  integral, 

oo 

(AT(f,l))  =  ^  J  dafvAA(a,t)  ,  (19) 

o 


where 

OO 

A(a,t)  =  j  pdpB(p,t)  . 

a 

The  time  for  diffusion  of  the  air  over  scales  of  the  order  of  a  is  tj)  =  a2/ D.  With  D  —  0.24 
cm2/sec,  a  =  0.2  ^m  one  has  <£>  =  1.7  nsec.  At  t  »  to  one  can  write 

A(a,t)  =  D  Ji  -  Ta  (l  -  .  (20) 

From  Eq.  (16)  the  thermal  blooming  absorption  coefficient  becomes^8) 

OO 

aa(F,t)  =  J  daf(a,F)aA(a )  (l  -  e_</r°)  (21) 

O 

By  using  the  distribution  function  (12)  one  can  find  aa  as  a  function  of  time  for  different 
LAV  laser  intensities  and  wavelengths.  The  results  are  plotted  in  Fig.  7,  where  the  absorp¬ 
tion  coefficient  is  normalized  to  that  in  the  absence  of  vaporization.  The  LAV  intensities 
axe  chosen  at  the  stimulated  rotational  Raman  threshold  for  atmospheric  propagation  to 
altitudes  of  10  km  (see  Section  2.6).  Note  that  at  short  times  a0  rises  linearly  in  time  and 
it  reaches  a  maximum  when  the  vaporization  threshold  is  reached.  After  that  it  falls  off 
exponentially  due  to  vaporization.  For  example,  a  KrF  laser  with  tp  —  1  psec  and  I  =  3 
MW/cm2  (fluence  of  3  J/cm2)  will  reduce  the  absorption  coefficient  by  a  factor  of  50. 
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NORMALIZED  ABSORPTION  COEFFICIENT,  «/« 


Figure  7:  Absorption  coefficient  in  air  due  to  aerosols 
17 
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2.5  PHASE  DISTORTIONS  DUE  TO  AEROSOL  VAPORIZATION 


The  heat  deposited  in  the  aerosols  and  transferred  to  the  air  changes  the  local  refractive 
index.  Let  us  define  S  =  8{n  —  l)/(n  —  1),  where  n  is  the  refractive  index  of  the  air.  S 
induced  by  the  LAV  pulse  is  given  by^2) 


1  d2  \  dS 
C2dt2)  dt 


3— i-VJ  [a.  (/,«)/] 
IP 


(22) 


where  Cs  is  the  sound  speed,  7  =  1.4  is  the  ratio  of  specific  heats,  p  is  the  ambient  pressure 
and  I  refers  to  the  intensity  of  the  LAV.  The  number  of  wavelengths  of  phase  distortion 
in  a  propagation  distance  L  is 

L 

8N  =  ^  f  Sdz  (23) 

J 

o 

Note  that  during  the  LAV  pulse  the  phase  distortion  8N  is  in  the  transient  thermal  bloom¬ 
ing  regime  and  it  is  negligible.  This  means  that  the  LAV  does  not  have  to  be  phase  com¬ 
pensated  for  thermal  blooming.  8N  is  carried  by  the  wind  from  the  LAV  guard-ring  into 
the  main  HEL  beam.  One  has  to  evaluate  8N  at  times  of  the  order  of  the  wind  clearing 
time  of  the  guard-ring,  tv  ~  dv/Vw.  tv  is  typically  0.1  sec  and  is  much  longer  than  the 
time  for  the  sound  wave  to  cross  the  guardring.  In  that  case  W2S  »  C~2d2S  and  the 
phase  distortion  becomes 


8N  =  -fc-  -  ,^7  --)■  [ dzl  dtaa(I,t)I 

A„7P  J  J 

o  o 

Prom  Eq.  (21)  aa  can  be  modeled  by  the  expression 

^  1  —  exp 


,r  .  (It 

aa(I,t)  =  a0exp[-y 


where  r  is  the  equilibration  time  during  vaporizations, 


(24) 


(25) 


r  =  raexp 


One  takes  into  account  the  exponential  decrease  in  the  absorption  coefficient  and  the 
effective  time  delay  (r)  for  the  transfer  of  heat  to  the  air.  As  seen  from  Fig.  7,  r  is 
typically  less  than  1  psec.  By  integrating  Eq.  (24)  with  aa  from  (25)  one  finds 


8N 


£ 

2ir 


£{i  _  (-^)  -  [l  -  exp  -  (L  +  i)  r,]  } 


where 


Ir  +  Fv 

27r(n  —  1)(7  —  1) 
Av7  P 


Q;o(A  »)P i? 
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When  tp  >  t,  It  <  Fv,  which  is  typical  for  the  vaporization  regime,  the  result  is  simplified, 

6N  =  ^L  1  -exp(-j£\  =  8N0  +  8Ni  (27) 

Significant  vaporization  occurs  when  Irp  >  Fv  and  the  last  term  in  Eq.  (27)  is  small. 
The  phase  change  is  nearly  constant  across  the  beam.  With  ot0  from  Eq.  (15),  3  x  103 
particles/cm3,  a0  =  0.1  /um  and  L  =  2  km  one  finds  the  LAV  induced  constant  phase 
distortion  at  different  wavelengths,  SN0  (10.6  pm)  =  0.15,  6N0  (1.06  /im)  =  1.5  and  SN0 
(0.248  p m)  =  6.  Clearly  the  KrF  laser  has  the  largest  phase  distortion,  but  even  6  waves 
are  negligible  for  the  HEL  phase  compensation  adaptive  optics.  A  different  problem  arises 
from  the  8Nj  term.  Even  when  the  spatially  averaged  SNi  «  SN0 ,  it  introduces  spatial 
scales  from  the  turbulence  induced  intensity  fluctuations  of  LAV.  These  scales  are  of  the 
order  of  the  coherence  length,  r0  ~  A5/6,  and  are  much  smaller  for  a  KrF  LAV.  The  actuator 
spacing  of  the  deformable  mirror  for  an  HEL  beam  at  A  =  1.06  pm  is  determined  by  r0 
which  is  typically  10  cm.  Such  a  mirror  does  not  compensate  at  scales  corresponding  to  a 

KrF  LAV  with  r0  =  3  cm  if  they  exceed  significantly  those  of  the  atmospheric  turbulence. 

However,  by  increasing  the  KrF  laser  fluence  the  deleterious  effects  of  small  scale  distortions 
can  be  minimized,  as  is  now  shown. 

The  figure  of  merit  is  the  Strehl  ratio  which  is  determined  by  the  mutual  coherence 
function, 

M12  =<  expi(ip(  1)  —  tp(2))  >  .  (28) 

The  phases  are  evaluated  at  two  points  separated  by  a  two-  dimensional  vector,  A.  <  ...  > 
means  ensemble  average.  The  turbulence  Mf2  is^11-1 

T  /A\5/3' 

Mf2  =  ezp  —3.44  f — J  ,  (29) 

where  the  coherence  length  is  r^3  =  0.423C3L(27r/A)2.  The  LAV  induced  mutual  coher¬ 
ence  function  from  Eq.  (27)  is 

L 

Mi 2  =<  J dz  | exp  —  exp  ^-I(r  +  A,z)^j  |  >  .  (30) 

O 

Note  that  in  deriving  Eq.  (27)  we  have  assumed  that  I  is  constant  in  z.  In  Eq.  (30)  we 
have  used  a  general  form  I(r,  z).  The  intensity  correlation  function  for  weak  turbulence  is 

(M(r, *,)«(?  +  A,  z,))  =  2.24CX'V‘/‘  ^1  -  0.6  (£^)  '  exp  (-^)  (I)2 

where 

1  27T 

81  =  I  -  (I),z  =  -(zi  +  z2)  ,£  =  zi  -  z2  and  kv  =  —  •  (31) 
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After  a  lengthy  calculation  one  finds  from  Eqs.  (30)  and  (31) 


When  the  fluence  is  such,  that  the  mutual  coherence  function  Mi  2  is  larger  than  the 
turbulence  M^,  the  LAV  induces  phase  changes  smaller  than  those  generated  by  turbulence 
alone.  In  that  case  the  HEL  compensation  system  will  work  and  the  LAV  phase  distortions 
can  be  neglected.  From  Eqs.  (29)  and  (32)  the  requirement  M12  >  Mj2  leads  to  the 
following  constraint 

0.2/SL^-exp  Jr)  <  1  •  (33) 

From  the  discussion  after  Eq.  (27)  note  that  (3L/ 2ir  =  6  at  \v  =  0.248  /xm.  At  A/A„ 
=  4  from  (33)  one  requires  that  F/Fv  >  5,  which  is  consistent  with  the  requirement  of 
a  significant  reduction  in  the  absorption  coefficient.  At  At,  =  1.06  /x m  the  corresponding 
requirement  would  be  F/Fv  >  1;  however,  at  A„  =  10.6  /xm  the  small  distortions  would 
not  be  an  issue. 

2.6  LASER  FOR  AEROSOL  VAPORIZATION  (LAV) 

In  this  section  we  describe  the  main  characteristics  of  a  laser  designated  to  clear  the 
atmospheric  path  of  a  HEL  beam  by  aerosol  vaporization.  Its  mission  is  to  reduce  the 
aerosol  absorption  by  at  least  an  order  of  magnitude  over  a  cross  section  of  a  few  square 
meters  and  up  to  altitudes  of  10-20  km.  This  will  permit  the  high  power  beam  for  BMD 
to  propagate  with  significantly  reduced  phase  distortion  due  to  thermal  blooming.  The 
intensity  of  LAV  is  limited  by  the  threshold  for  stimulated  rotational  Raman  scattering  in 
N2.  The  SRRS  steady  state  gain  coefficient  is  given  approximately  by  g  (cm/TW)  =  2.5/A 
(jum).  Collisional  broadening  dominates  over  Doppler  broadening  up  to  altitudes  of  50  km; 
therefore,  to  a  good  approximation  g  may  be  considered  to  be  constant  at  altitudes  below 
50  km,  and  zero  beyond.  For  an  integrated  gain,  G  =  glthh  —  30  from  the  ground  up  to  an 
altitude  of  h  =  10  km,  the  threshold  intensity  is  Ith  (MW/cm2)  =  12  A  (/xm).  In  the  case 
of  a  CO2  wavelength  one  might  not  be  able  to  operate  near  the  Raman  threshold  because 
of  optical  breakdown  in  air  (due  to  aerosols).  The  fluence  of  the  LAV  has  to  exceed  Fv 
to  vaporize  the  particles,  and  it  has  to  deliver  the  energy  for  pulselengths  tp  <  r0,  the 
particle  equilibration  time.  This  will  insure  effective  heating  of  the  aerosols  and  not  the 
surrounding  air.  The  requirement  for  the  LAV  pulselength  is  ra  >  rp  >  tc  =  Fv/Ith- 
Note  that  ra  =  1.3  /xsec  for  the  dominant  particles  with  radius  a  —  0.2  /xm.  At  the  KrF 
wavelength  Fv  =  0.8  J/cm2  and  rc  =  0.3  /x sec.  For  A  =  1.06  //m,  Fv  =  12  J/cm2  and  rc 
=  1  /xsec.  At  the  CO2  wavelength  Fv  =  8  J/cm2  and  rc  =  7  x  10~2  /xsec.  Even  at  J(10.6 
/xm)  =  40  MW/cm2  one  has  rc  =  0.2  /xsec,  which  may  be  a  reasonable  regime  to  operate 
to  avoid  optical  breakdown.  Clearly,  the  KrF  and  CO2  wavelengths  are  well  suited  for 
aerosol  vaporization,  while  the  1  /xm  wavelength  is  a  borderline  case. 

As  a  spatial  format  for  LAV,  consider  a  guard-ring  of  thickness  d„,  surrounding  the 
main  HEL  beam  of  diameter  d  (Fig.  8).  In  this  format  the  beam  will  be  cleared  indepen¬ 
dently  of  the  wind  direction  and  the  wind  shear.  The  area  cleared  is  ( d  »  dv)  Av  ~  7r ddv. 
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The  pulse  repetition  rate  of  the  LAV  is  determined  by  the  wind  velocity  and  dv,  PRF  = 
Vw/dv.  For  V\y  —  5  m/sec  and  d v  =  0.5  m  one  finds  PRF  =  10.  At  this  repetition  rate  all 
aerosols  moving  through  the  guard-ring  axe  irradiated  by  at  least  one  LAV  pulse.  Thus, 
after  a  wind  clearing  time,  the  whole  FEL  beam  will  propagate  at  a  reduced  absorption  in 
steady  state.  The  average  power  of  the  LAV  can  be  estimated  as 

Pv  —  it  dr pVwI  ,  (34) 

where  I  is  the  peak  LAV  intensity.  Since  Irp  >  Fv,  for  vaporization  one  finds 

Pv>7cdVwFv  (35) 

For  an  FEL-TIE  beam  diameter  of  4  m  one  finds  Pv  >  0.5  MW  at  A  =  0.248  pm.  At  A 
=  1.06  pm,  Pv  >  7.5  MW  and  for  A  =  10.6  /tm,  Pv  >  5  MW.  Again  the  KrF  wavelength 
is  the  most  promising,  but  even  at  the  other  wavelengths  the  power  of  the  LAV  is  a  small 
fraction  of  the  100  MW  plus  GBL  power.  The  scheme  of  using  LAV  for  clearing  up  the 
beam  path  appears  to  be  very  cost  effective. 

2.7  SCALING  OF  THE  CRITICAL  HEL  POWER  WITH  LAV 


The  atmospheric  absorption  as 


a  function  of  altitude  is  of  the  form 


a  —  a 


5 


(36) 


where  L  is  of  the  same  order  as  the  thickness  of  the  turbulence  layer,  Lt  —  L  —  2  km. 
The  strength  of  the  thermal  blooming  is  described  by  a  distortion  number, 


Nd  —  47r 


(n  —  1)(7  —  l)a(o)X 
•ypVw^d 


(37) 


where  P  is  the  average  laser  power,  d  is  the  beam  aperture,  n  is  the  refractive  index,  n  — 1  = 
3  x  10-4  at  1  atm  and  p  is  the  ambient  pressure.  For  laser  powers  in  the  MW  range  and 
beam  diameters  of  a  few  meters  Njj  is  typically  Nd  >  102.  The  phase  compensation  system 
can  give  Strehl  ratios  approaching  one  only  when  Nd  is  below  a  certain  threshold  value, 
N Dth .  When  one  neglects  wind  shear,  Nnth  can  be  evaluated  numerically  and  analytically, 
taking  into  account  phase  compensation  for  atmospheric  turbulence  and  thermal  blooming. 
The  result  is 

Nd  ~  1 61V|/12  ,  (38) 

where  Nt  =  r2/A L  is  the  turbulence  coherence  length  Fresnel  number.  The  threshold 
distortion  number  from  Eqs.  (37)  and  (38)  leads  to  a  critical  laser  power  for  the  propagation 
of  a  ground  based  laser, 


7pIVA7/127'o^6d 

(n  —  1)(7  —  l)L17/12a(o) 


(39) 


At  a  given  aperture  and  wavelength  only  laser  powers  P  <  Pc  can  propagate  through  the 
atmosphere  without  serious  degradation  in  optical  beam  quality.  Note  that  Pc  oc  d,  a  result 
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TEMPORAL  FORMAT  OF  LAV  AND  FEL 

FEL  PULSES  @  15-30  kHz 


SPATIAL  FORMAT 


LAV  BEAM 
FEL  BEAM 


Figure  8:  Temporal  and  spatial  pulse  format  for  vaporizing  aerosols 
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which  is  very  unfavorable  for  scaling  to  the  large  powers  necessary  for  BMD  applications 
(P  >  100  MW).  At  A  =  1.06  /zm,  r0  =  10  cm,  Vw  =  5  m/sec,  L  =  2  km,  a(o )  =  2x  10"7 
cm-1,  the  critical  power  is  plotted  as  a  function  of  beam  diameter  in  Fig.  9.  Note  that  at 
d  —  10  m  only  10  MW  of  power  can  be  transmitted  through  the  atmosphere  with  Strehl 
ratios  close  to  1.  More  realistic  models,  including  the  effect  of  wind  shear  appear  to  give 
an  increase  of  Pc.  The  scaling  in  that  case  is  still  the  subject  of  an  active  study.  On  the 
other  hand,  a  defensive  system  should  be  robust  and  work  well  under  stressful  conditions 
and  it  should  not  rely  on  favorable  weather.  We  consider  the  scaling  (39)  a  conservative 
estimate  which  one  might  be  forced  to  accept  as  a  guide  for  the  design  of  a  GBL. 

The  only  certain  way  to  increase  the  critical  power  is  to  reduce  significantly  the  absorp¬ 
tion  coefficient.  This  can  be  achieved  with  a  LAV  of  modest  powers.  Note  the  exponential 
dependence  of  a  on  LAV  fluence,  Eq.  (14).  By  using  the  power  of  LAV  in  a  guard-ring, 
Eq.  (34),  one  finds  the  aerosol  absorption  coefficient  with  vaporization. 


Now  one  can  replace  cv(o)  in  Eq.  (39)  with  the  reduced  value  a,.  The  scaling  law  of  the 
critical  HEL  power  becomes  a  function  of  beam  diameter  and  LAV  power.  The  result  for  a 
KrF  LAV  is  shown  on  Fig.  9  at  Pv  =  2  MW.  Note  the  enormous  (>  100  MW)  HEL  powers 
at  small  beam  diameters.  Unfortunately,  in  practice,  even  in  the  presence  of  LAV  one  can 
not  propagate  such  large  laser  powers  through  the  atmosphere  because  of  the  inherently 
low  duty  factors  of  the  HEL  and  the  threshold  for  stimulated  rotational  Raman  scattering 
(SRRS)  mentioned  in  Section  2.6.  At  duty  factors  <  10-3,  1  kW/cm2  of  average  intensity 
leads  to  peak  intensities  in  excess  of  1  MW/cm2.  The  HEL  is  constrained  by  a  Raman 
threshold  for  the  atmosphere  to  altitudes  of  50  km,  not  the  10  km  altitude  sufficient  for 
the  LAV.  The  peak  intensity  threshold  for  an  HEL  is 

J„<UW/«m*)  =  jgL  .  (41) 

In  terms  of  average  power  the  threshold  becomes 

Pth(UW)  =  -^Td2(m)rL(/zsec)PRF(kHz)  .  (42) 

A(jumj 

The  units  for  the  various  quantities  are  noted  in  brackets.  For  a  given  beam  diameter 
only  P  <  Pth  can  propagate  without  degradation  from  Raman  scattering.  A  HEL  with  tl 
=  50  nsec,  PRF  =  5  kHz  leads  to  the  Raman  threshold  curve  shown  on  Fig.  9  with  the 

“forbidden”  region  shaded.  By  way  of  example,  a  KrF  LAV  at  2  MW  of  power  can  clear 

a  channel  for  the  propagation  of  100  MW  of  HEL  power  at  beam  diameters  of  4-5  m. 
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3.0  SBIR  PHASE  II  PROGRAM 


SRL  will  propose  an  SBIR  Phase  II  program  to  determine  the  reduction  of  the  atmo¬ 
spheric  aerosol  absorption  by  aerosol  vaporization.  The  success  of  this  program  will  have  a 
significant  and  immediate  impact  on  the  design  of  a  ground  based  laser  for  ballistic  missile 
defense. 
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CONCLUSION 


As  has  been  discussed  in  this  report,  propagation  of  ground-based,  high-energy  laser 
(HEL)  beams  through  the  atmosphere  at  a  laser  wavelength  of  A i  =  1.06  //m  is  limited  by 
thermal  blooming  due  to  aerosol  absorption.  The  scaling  of  critical  laser  power  with  beam 
diameter  requires  that  the  aperture  size  exceed  20  m  for  average  laser  powers  of  100  MW. 

In  this  SBIR  Phase  I  report,  Science  Research  Laboratory  has  shown  that  this  problem 
may  be  solved  by  vaporizing  the  aerosols  with  a  designated  laser  for  aerosol  vaporization 
(LAV).  The  dependence  of  vaporization  on  fluence  and  wavelength  was  determined.  In 
particular,  it  was  shown  that  a  LAV  wavelength  \v  =  1.06  //m  and  intensity  of  10  MW/cm2 
yields  a  reduction  by  a  factor  of  5  at  2  //sec,  corresponding  to  20  J/ cm2  of  laser  fluence  per 
pulse.  The  best  LAV  wavelength  appears  to  be  that  of  a  KrF  laser,  A„  =  0.248  //m.  For 
this  wavelength,  an  intensity  of  3  MW/cm2  and  a  pulse  duration  of  1  //sec,  corresponding 
to  a  fluence  of  3  J /cm2  per  pulse,  reduces  aerosol  absorption  by  a  factor  of  50.  Such  a 
decrease  in  the  absorption  coefficient  eliminates  the  difficulties  associated  with  thermal 
blooming. 

The  dependence  of  the  critical  HEL  power  on  the  beam  aperture  diameter  d  has  been 
determined  in  the  presence  and  absence  of  LAV,  assuming  no  wind  shear.  Without  vapor¬ 
ization  the  scaling  law  is  pessimistic — one  can  propagate  less  than  10  MW  average  power 
at  d  =  10  m.  If  one  uses  a  KrF  laser  at  an  average  power  of  2  MW,  a  pulsewidth  of  1 
//sec,  and  a  pulse  repetition  frequency  (PRF)  of  10  Hz,  the  critical  power  at  d  —  5  m  is 
100  MW.  This  greatly  alleviates  the  thermal-blooming  problem. 

Science  Research  Laboratory  will  shortly  submit  a  Phase  II  SBIR  proposal  to  experi¬ 
mentally  test  and  demonstrate  the  LAV  concept. 
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